The first bifurcation point for Delaunay nodoids 

Wayne Rossman 

Abstract: We give two numerical methods for computing the first bifurcation point for 
Delaunay nodoids. With regard to methods for constructing constant mean curvature 
surfaces, we conclude that the bifurcation point in the analytic method of Mazzeo- 
Pacard is the same as a limiting point encountered in the integrable systems method of 
Dorfmeister-Pedit-Wu. 

1. Introduction 

Delaunay surfaces in Euclidean 3-space 

R 3 = {(xi,x 2 ,x 3 ) | Xj € R} 

are constant mean curvature (CMC) surfaces of revolution, and they are translationally periodic. By a 
rigid motion and homothety of R we may place the Delaunay surfaces so that their axis of revolution 
is the xi-axis and their constant mean curvature is H — 1 (henceforth we assume this). 

We consider the profile curve in the half-plane {(xi, 0, £3) € R 3 | X3 > 0} that gets rotated about 
the xi-axis to trace out a Delaunay surface. This curve alternates periodically between maximal and 
minimal heights (with respect to the positive X3 direction), which we refer to as the bulge radius and 
the neck radius, respectively, of the Delaunay surface. Let us denote the neck radius by r. 

Delaunay surfaces come in two 1-parameter families: one is a family of embedded surfaces called 
unduloids that can be parametrized by the neck radius r G (0, 1/2]; the other is a family of nonembedded 
surfaces called nodoids that can be parametrized by the neck radius r £ (0, 00). For unduloids, r = 1/2 
gives the round cylinder. For both unduloids and nodoids, the limiting singular surface as r — > is a 
chain of tangent spheres of radii 1 centered along the xi-axis. 

In this paper we shall be concerned with nodoids. We will see that a common bifurcation point for 
Delaunay nodoids is encountered in the following two distinctly different approaches for constructing 
CMC surfaces: 

(1) Using analytic techniques, Mazzeo and Pacard [15] showed existence of a finite value tq so that 
for neck radii r < ro the nodoids are nonbifurcating, and for r > ro the nodoids can bifurcate. 
Bifurcating nodoids are of interest because they deform smoothly through families of CMC 
surfaces that are of bounded distance from a fixed line yet are not surfaces of revolution. Before 
the work [15], such examples were unknown. 

To study this bifurcation, a particular Jacobi operator associated to the second variation 
formula for Delaunay surfaces is used, along with a function space that contains only functions 
with the same translational periodicity as the Delaunay surfaces themselves. Because these 
functions are translationally periodic, they do not have finite L 2 norms on the entirety of 
the surfaces. Hence bifurcation is a different notion from that of " nondegeneracy" of CMC 
surfaces, i.e. CMC surfaces with no nonzero Jacobi fields of finite L 2 norm on the entire 
surfaces. However, the two notions have the common trait of being highly useful tools for 
producing previously unknown CMC surfaces. (There has been much work done related to the 
nondegeneracy of CMC surfaces and construction of new CMC examples, see for example the 
works of Kapouleas, Kusner, Mazzeo, Pacard, Pollack, Ratzkin [10], [11], [14], [16], [17], [21].) 
As our interest here is in the bifurcation of Delaunay nodoids, from the outset we consider only 
periodic functions on Delaunay surfaces. 

Mazzeo and Pacard gave a clear reason for the existence of this bifurcation point ro, in terms 
of the existence of nontrivial nullity for the Jacobi operator, but they did not compute the 
precise value of ro- 

(2) Using integrable systems techniques developed by Dorfmeister, Pedit and Wu in [8] , Dorfmeister, 
Wu [9] and Schmitt [23] (see also [12]) constructed genus CMC surfaces with three asymp- 
totically Delaunay ends. In [9] the construction was restricted to surfaces with asymptotically 
unduloidal ends, because such ends are embedded. However, the construction in [23] and [12] 
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includes asymptotically nodoidal ends as well. The construction begins with the selection of 
a certain "DPW potential", and the DPW potential in [23] and [12] fails to exist when some 
nodoidal end has an asymptotic neck radius greater than 1/2. Furthermore, the arguments in 
[23] and [12] showing that each end converges to a Delaunay surface work only when the limiting 
Delaunay surface is an unduloid, or is a nodoid with neck size less than 1/2. This suggests that 
there is possibly some obstruction in the DPW approach to CMC surfaces with Delaunay ends 
that occurs only for asymptotically nodoidal ends with asymptotic neck radii at least 1/2. 

In the construction in [23] and [12], this limiting value 1/2 appears explicitly, but the under- 
lying reasons for its appearance are left unexplained. 

It is natural to ask if the bifurcation radius r$ in the first approach coincides with the limiting value 
1/2 in the second approach. The purpose of this article is to numerically confirm this: 

Numerical result. The bifurcation radius r is equal to 1/2. 

This result is of interested with respect to both approaches above. For the first approach, it gives 
the exact (previously unknown) value for tq. For the second approach, it provides a reason (via the 
bifurcation properties shown by Mazzeo and Pacard) for the existence of the previously mysterious 
limiting value 1/2. 

In fact, we shall show a stronger numerical result about the first eigenvalue of a particular operator, 
for which the above numerical result is an immediate corollary. 

To provide added confidence in the accuracy of our numerical arguments, we give two different 
independent algorithms for showing that r = 1/2 (and for showing the stronger numerical result as 
well). The first method involves using an ordinary differential equation (ODE) solver and symmetry 
properties of the first eigenfunction. The second method requires more machinery, the basic tool being 
the Rayleigh quotient characterization for eigenvalues of a self-adjoint operator. It involves numerical 
integration of smooth bounded functions of a single variable on a finite interval, and it has the advantage 
of also giving estimates for other eigenvalues beyond the first one. 



2. Jacobi elliptic functions 

Before introducing equations for nodoids in the next section, we first briefly review the properties of 
Jacobi elliptic functions needed here. For (j),k £ R with < k 2 < 1, we set 

r<t> 

a) a 



o y/l — k 2 sin 2 tp 



and then we define 



sn k (£)=sm(fi, cn fc (£) = cos</> , dn fe (£) = \J I - k 2 si 



sin 2 i 



Extending these functions analytically to £ e C, they are defined on the complex plane (with singular- 
ities). A short computation gives 

(2) ^sn fc (£) = cn fc (£).dn fc (£). 
Defining 

(3) z = c ■ sn fc (a£ + b) 

for constants a, b and c, Equation (2) then gives 

, . 1 dz _ a ■ cn fe (a£ + b) ■ dn fc (a£ + b) 

zdt; sn fe (a£ + 6) 

and 

(5) (^j 2 + a 2 (i + kfz 2 - a 2 (5 + k~ c - x z 2 f = o . 



Let sn fe 1 denote the (multi- valued) inverse function of sn k . Although multi- valued, from here on out 
let us fix sn^T 1 (l) to be equal to the specific value of £ e R given by the integral (1) with (f> = ir/2. 

Lemma 1. Suppose that b =sn^ 1 (l), and that c > and a € iR. Then the function z as in Equation 
(3) is positive (i.e. z G R + ) and periodic with respect to the variable £ G R. 

Proof. Using the identity 

sn k (u)cn k (v)dn k (v) + sn k (v)cn k (u)dn k (u) 



sn fe (u + v) 



1 - k 2 sn 2 (u)sn 2 (v) 
with u = a£ and v = b, we have 

.cn fc (a£) 
z = c- 



~dn fe (a£) " 
Then using the relations 

f~\ 1 A dn fc'(^) 

cn k (iu) = -— , dn k {iu) = — 

cn fc / (u) cn fc / (u) 

with k' e R satisfying (fc') 2 + k 2 = 1, we have 



z = 



dn fe /(-m£) 

It follows that z is positive and periodic when £ S R, because — ia£ € ii and dn^/ is positive and 
periodic with respect to a real variable. □ 

3. Parametrizing nodoids 

Let (xi, X2, X3) be the usual rectangular coordinates for R 3 . Consider a Delaunay nodoid with the 
xi-axis as its axis and with constant mean curvature H = 1. Let 

(x(t),z(t)), teR 

be a parametrization of the profile curve of the nodoid in the xix 3 -plane, and so the surface can now 
be parametrized by 

V(t,6) = (x{t),z{t) cos 6, z(t)sm6) , teR, 6 e [0,2tt) . 

Suppose further that the parameter t is chosen to make the mapping T>(t, 9) conformal with respect to 
the coordinates (t, 9). Let t = a and t — b be values at which the nodoid achieves two adjacent necks, 
i.e. z(t) has local minima at both t — a and t = b equal to the neck radius r but at no t € (a,b). 
Conformality implies that the first fundamental form is 

ds 2 = {{x') 2 + (z') 2 )dt 2 + z 2 d9 2 = p 2 (dt 2 + d9 2 ) , 

with 

p 2 = {x' f + (z 1 ) 2 = z 2 . 

The second fundamental form is then 

\{x"z' -z"x')dt 2 + x'd6 2 , 

and so the coordinates (t, 9) are curvature line coordinates, that is, the coordinates are isothcrmic. 
Furthermore, the mean curvature H = 1 implies 

2z 3 - z'x" + x'z" - zx' = . 

This has a first integral, that is, using z 2 = (x 1 ) 2 + (z') 2 and zz' = x'x" + z'z", it is equivalent to 

{Ax' - Az 2 )' = . 

Thus 

(6) to = Ax' - Az 2 
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is satisfied for some i-indcpcndent constant m. Because there are points on a nodoid where x' = 0, we 
have 

m < . 

In fact, m e (— oo, 0) is a parameter that determines the full family of Delaunay nodoids. Equation (6) 
can also be established using a homology invariant on CMC surfaces called the weight, as explained by 
Korevaar, Kusner and Solomon in [13]. We define m as the mass: 

Definition 1. Given a nodoid T>(t, 9) parametrized as above satisfying Equation (6), we say that m is 
the mass (also sometimes called the weight or flux) of the nodoid T>(t, 9). 

Equation (6) implies 
(7) z' = ±^z 2 - (f + z 2 ) 2 , m<0. 

Then x' is determined by z via Equation (6), so x is determined up to translation, as expected. 
It follows from Equation (5) that Equation (7) has the solution 

z = -2B ■ sn B/A ( X (t)) , X (t) = 2iA(t - b) + sn^l) , 

where B = — ni — 1) and A = ^ — B. Lemma 1 shows us that this solution is positive and 

periodic with respect to i, so it is indeed the height function for the profile curve of a nodoid in the 
upper half of the xiX3-plane. 

From this we can see that the minimum (neck radius) of z is 



(8) 



rn 



2 

at t = a and t — b, the maximum (bulge radius) of z is 



&tt = z±k, and x' = and z = yj =f at both t = ^ and t = 
The Gaussian curvature K is determined by ds 2 as 

P 2 

where 

(9) A = K + K 

is the standard Euclidean Laplacian operator. 

We now introduce a function V that will be used in the second variational formula later. We set 

V = (4 - 2K)p 2 . 

It can be computed that 

(10) V = z- 2 {2z i + \m 2 ) . 
Note that V = V(t) is a function of t, and is independent of 9. 
Lemma 2. The function V = V(t) has the following symmetries: 

V(£ + t) = V(t-t) Vie {a, ^,^,^,6} . 

Furthermore, 

-m < V < 2 - m MteR. 
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FIGURE 1. A nodoid and a CMC surface with three ends asymptotic to nodoids. One of two 
halves of each surface is shown here. (The computer graphics were made using N. Schmitt's 
cmclab software [24] and K. Polthier's Javaview software [20].) 

Proof. Wc will show that V(a + t) = V(a-t) and V(^d± + t) = V('^ - t), and then all other 
symmetries in the lemma follow. Because z(t) itself has the symmetry z(a + t) = z(a — t), clearly also 
V(a + t) = V{a-t) by Equation (10). To show V(^£ + t) = V(^^--t), we first note, using Equation 
(4), that 

z'(t) = 2iA ™B/A (x(t))dn B/A ( X (t)) 
z(t) sn B/A ( x (t)) 

and so 

Substituting m 2 /8 = (2 - m)z 2 - 2z 4 - 2(z') 2 from Equation (7) into Equation (10), wc have 
(12) V(t) = 2-m-2^y . 

Then by Equation (11), we have V(^^- +t) = V{^^- - 1). Equation (12) also implies that V < 2 - m 
for all t. 

To see that V > — m for all t, we simply consider V = 2( + as a function of ( = z 2 > 0. 

2 

Elementary calculus gives that 2£ + ^ > — m for all Q > 0. □ 



4. Second variation for nodoids 

We now consider an arbitrary volume-preserving periodic smooth variation of the surface, that is, a 
family of surfaces T> s (t, 9) that is 

(1) smooth with respect to s € (— e, e) for some e > 0, 

(2) an immersion with respect to the coordinates (t, 9) for any fixed s € (— e, e), 

(3) satisfying V (t, 9) = V(t, 9) for all f € 12 and 9 G [0, 2tt), 

(4) periodic with the same period for all s £ (— e, e), i.e. T> s (t%, 9) = T> s (t2, 9) for all s S (— e, e) and 
all t% — t\ an integer multiple of b — a, 

(5) with compact regions bounded by the surfaces T> s (t, 9),t G [a, b], 6 G [0, 27r) and the two disks 

D a (s) = {(x(a),x 2 ,x 3 ) \x 2 2 +x 2 3 < (z(a) + O(s)) 2 } , 

D b (s) = {(x(b),x 2 ,x 3 ) \x 2 2 +xl< (z(b) + O(s)) 2 } 
having the same volume for all s G (— e, e). 
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The disks D a (s) and D (s) are allowed to deform smoothly in s, and are not necessarily perfectly 
round when s ^ 0. Thus these deformations do not have any fixed "boundary" curves. The essential 
properties of these deformations are only that they are both periodic and volume-preserving. 
Let TV denote a unit normal vector to T>o(t, 9). We define the function u by 

u = u(t,0) = (£v s (t,6)\ s=0 ,N)eF. 

The volume-preserving condition 5. above implies 

r b r 2ir 

up 2 d9dt = . 



// 

Ja JO 



Definition 2. We will call the compact portion T>(t, 0), t £ [a, b], 6 £ [0, 2tt) a fundamental piece of the 
nodoid. 

Let A(s) denote the area of the surface T> s (t, 9) for t £ [a, b], 9 £ [0, 2ir). Because the nodoid is CMC 
H = 1, the first variation formula for a fundamental piece is 



d At \ 



-fl 

=0 J a JO 



2-rr 

2 



up A d9dt = 



Thus it is the second variation formula for volume- preserving variations (see [1]) 

d2 At , 
d^ A{3) 



nO nZTZ 

= jj u(-A ds 2u- (4-2K)u)p 2 d9dt 

J a JO 

nb nl-K 

\ / u ■ C{u)d9dt = , C(u) := -Au-Vu , 

Ja JO 



with A ds 2 (respectively A) the Laplace-Beltrami operator determined by ds 2 (respectively the Euclidean 
Laplacian as in (9)), that will determine if the variation increases or decreases area (when this formula 
is nonzero). We are applying the first and second variation formulas only to periodic variations here, 
but we remark that they hold for other nonperiodic types of variations as well. 

5. Spherical harmonics 
We now consider the eigenvalue problem for C. We first introduce the function space 
T = {u = u(t, 9) e C°°(t, 9) | u(t u 9) = u(t 2 ,9) for h - t 2 e {b - a) ■ Z}. 
The eigenvalue problem is to find u £ T and A € R such that 
(13) C(u) = Xu . 

Let us decompose such an eigenfunction u into its spherical harmonics: u = u(t, 9) can be written as 

u = u (t) + ^2(u ji+ (t) ■ cos(jff) + uj-(t) ■ sin(j0)) , 

where uo(t), Uj.+(t) and Uj t -(t) are periodic functions of t, that is, they lie in the smaller function space 

T = ?T\{u = u(t, 9)1^ = 0} , 

i.e. those functions in T that are independent of 9. Defining the operators 

C - - V + i 2 
'-j — at 2 v + J 

for j e Z + U {0} on the function space P, the relation (13) gives also that 

^o(uo) = Xu , Cj{uj t ±) = Xuj t ± 

for all j > 1. So the following lemma is immediate: 

Lemma 3. A real number A is an eigenvalue of C, i.e. C(u) = Xu for some u £ T , if and only if X is 
also an eigenvalue of Cj for at least one value of j £ Z + U {0} 7 i.e. Cj{u) = Xu for some u £ T and 
some j £ Z + U {0}. 
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The following lemma is analogous to Proposition 4.4 in [15], but different notations were used there: 
Proposition 4. Both —1 and are eigenvalues of Co- 

Proof. We first show that is an eigenvalue. Considering the vector (1,0, 0) as a constant vector field of 
R 3 , the associated translational flow gives a periodic volume-preserving deformation T> 8 of the nodoid 
V = Vq. Because this flow is actually a family of rigid motions, it is also area-preserving. It is a 
classical fact that the scalar product of the unit normal vector of a CMC surface with a Killing field is 
in the null-space of the Jacobi operator of the surface (see [7], page 196, or the proof of Theorem 2.7 in 
[2]). Hence, if (1, 0, 0) is decomposed into (1, 0, 0) = uN + v at each point of V with uN normal to T> 
and v tangent to V, and with u € T and \N\ = 1, then C(u) = 0. Since u € T , it follows that is an 
eigenvalue of Co for any choice of m < 0. Also, by direct computation, we have u = —z'/z, and then 
using (z') 2 = z 2 — (z 2 + m/4) 2 and z" = z — 2z(z 2 + to/4), we have Cq(u) = 0. 

Now we show that —1 is an eigenvalue. Similarly to the previous case with eigenvalue 0, we now 
consider the vector (0, 0, 1), producing a smooth vector field of B? , whose associated translational flow 
again gives a periodic volume-preserving deformation T> s of T> = T>o that is once again area-preserving. 
Now we decompose (0,0,1) = u sin ON + v at each point of V with u sin ON normal to V and v tangent 
to V. Then C(u sin 9) = 0, and u e T is an eigenfunction of C\ with eigenvalue 0. It follows that —1 is 
an eigenvalue of Co for any choice of m < 0. Again, we could also see this by direct computation: we 
have u — x'/z = z + to/(4z), and then using (z') 2 = z 2 — (z 2 + to/4) 2 and z" = z — 2z(z 2 + m/4), we 
have Cq(u) = —u. □ 

Both the operator C on the function space T and the operators £j on the function space T are 
essentially self-adjoint, and hence standard functional analysis arguments (see [3] or [26] or [5] for 
example) give the following result: 

Proposition 5. Each of the operators L on the function space T and Cj on the function space T 
satisfy the following: 

(1) all eigenvalues are real, and there are a countably infinite number of them, 

(2) all eigenvalues are greater than some real constant, 

(3) the eigenvalues do not accumulate at any finite real value, 

(4) when written in increasing order, the eigenvalues increase to +oo, 

(5) the eigenspace associated to each eigenvalue is finite dimensional, 

(6) the collection of eigenspaces spans the full function space. 

We refer to the eigenvalue that is less than all the others as the first eigenvalue. The second eigenvalue 
is the one that is less than all but the first eigenvalue. The third eigenvalue is the one that is less than 
all but the first and second eigenvalues, and so on. 

The Courant nodal domain theorem, which is valid in our setting (see [6]), tells us that the number of 
nodal domains of the eigenfunction associated to the fc'th eigenvalue is at most k (here we are counting 
eigenvalues with multiplicity). If the eigenspace associated to the first eigenvalue contained two linearly 
independent functions, then some linear combination of the two functions would be a function that 
attains both positive and negative values. This would contradict the Courant nodal domain theorem, 
so the eigenspace associated to the first eigenvalue is 1 dimensional. 

Also, in the case of the operator Cj, the domain on which the functions are defined is 1 dimensional, 
so Cj is an ordinary differential equation, not a partial differential equation. Hence, saying that any 
(not identically zero) eigenfunction associated to the first eigenvalue has at most one nodal domain is 
equivalent to saying that any such eigenfunction never attains the value zero. We conclude the following: 

Proposition 6. For both the operator C defined on the function space T and the operator Cj defined on 
the function space T , the eigenspace associated to the first eigenvalue is 1 dimensional. Furthermore, 
in the case of Cj, any not-identically-zero function in this eigenspace never attains the value zero. 

Remark. The eigenvalues —1 and of Co are actually the second and third eigenvalues of Co- This was 
shown in [15] by counting the nodal domains of the corresponding eigenfunctions. 
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FIGURE 2. Plots of the functions /i = m, fz = m — 2 and /2 = Vdt. Ao is also a 

function of ra, and Lemmas 7 and 9 imply that Ao must lie between f^ and fa. 



6. Primary result 

The first two operators £o an d C\ never give any bifurcation in the sense of [15]. The first bifurcation 
point (as in [15]) is defined as follows: 

Definition 3. The first bifurcation point for Delaunay nodoids occurs at the value of m < closest to 
zero for which Cj for some j > 2 has a nonpositive eigenvalue, i.e. 

CjU — Xu with A < and u e T for some j > 2 . 

This occurs at the largest value of to for which the first eigenvalue of £2 is zero, and is equivalent 
to the first eigenvalue A of £ being —4. This A can be computed using a minimum of the Rayleigh 
quotient ([3], [26], [5]): 



(14) A = 



mm 



J b a u(t)C (u(t))dt 



u(t)er\{o} f*u{t) 2 dt 



Remark. Bifurcation points are points where bifurcation actually occurs, that is, in any neighborhood 
of the Delaunay surface there are CMC surfaces that are not rotationally symmetric. What we have 
defined in Definition 3 is an "infinitesimal symmetry breaking point" where bifurcation might happen, 
because having zero as an eigenvalue of £2 is a necessary condition for the bifurcation to occur. However, 
because the multiplicity of the zero eigenvalue of £2 is one (see Proposition 6), and because the derivative 
of the first eigenvalue with respect to the parameter to is not zero (see the primary numerical result 
below), it turns out that the infinitesimal symmetry breaking point is an actual symmetry breaking 
point, i.e. a true bifurcation point [19]. 

We will give our primary numerical result and two different numerical methods for showing it. But 
first let us give two simple mathematically rigorous lemmas that support this upcoming numerical 
result. Both of these results are immediate from the Rayleigh quotient formulation for Ao above. The 
first lemma is also shown in Proposition 4.4 of [15], with different notations. 

Lemma 7. For a nodoid of mean curvature 1 and mass m, the first eigenvalue Ao of Cq satisfies 
m — 2 < A < m. 

Proof. Using the Rayleigh quotient characterization (14), this follows directly from the fact that —to < 
V = V(t) < 2 — m for all t, as shown in Lemma 2. □ 

We then immediately have: 

Corollary 8. For nodoids of mean curvature 1 and mass to, the first bifurcation point occurs for some 
to between —4 and —2. 
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Lemma 9. For a nodoid of mean curvature 1 and mass m, the first eigenvalue Ao of Co satisfies 

Ao<^r/>dt. 

Proof. Inserting u = 1 into the Rayleigh quotient in Equation (14), this lemma follows. □ 

Numerically checking that J a Vdt is an increasing function of m < that becomes —4 when (and 
only when) to is approximately —3.036, we have the following: 

Preliminary numerical result. For nodoids of mean curvature 1 and mass m, the 
first bifurcation point occurs for some m > —3.036. 

Now, armed with the knowledge that the bifurcation point must occur for some m between —3.036 
and —2, we state our primary numerical result: 

Primary numerical result. For a nodoid of mean curvature 1 and mass to, the first 
eigenvalue Ao of £ is 

Ao = to — 1 . 

Hence the first bifurcation point for nodoids occurs when m = — 3. 

By the formula (8), m = — 3 precisely when the neck radius is r = 1/2, thus the numerical result 
stated in the introduction is a direct corollary of this numerical result here. 

7. First method 

In this and the next section we give two independent methods for numerically confirming our primary 
numerical result. The first method in this section is the simpler of the two. 
We wish to solve 

d 2 

(15) -— u + (V + X)u = 

for some A e [m — 2, to] and some function u = u{t) € T . If there is only one such value for A < — 1 in 
the range [to — 2, m] (and we will see that this is the case), then this A will be the first eigenvalue Ao- 
We begin with a mathematically rigorous lemma that this first numerical method is based on: 

Lemma 10. We fix any m < and consider the eigenvalue equation (15). Let u be a not-identically- 
zero eigenfunction associated to the first eigenvalue Ao . Then u cannot attain zero at any value of t, 
and 

u(£ + t) = u(£-t) for all £ E {a, ^±k, 2±5*, b} . 

In particular, 

u'(a) = u'(^) = u'(2±*) = u'(2±2*) = u'(b) = . 

Proof. Proposition 6 tells us that the eigenspace associated to the first eigenvalue Ao must be one 
dimensional. If u did not have all of the same symmetries as V as in Lemma 2, then there would be 
some t € {a, ^j^-, ^rp, ^^,6} so that u(£ + t) and u(£ — t) are two linearly independent eigenfunctions 
associated to the same first eigenvalue Ao- This is a contradiction, implying the lemma. □ 

Because multiplying u by a real scalar does not affect Equation (15), we may assume that 

= i. 

The numerical method is then to numerically solve Equation (15) for u with initial conditions 

U (32±6 ) = lj u '(32±*) =0 

by some ODE solver (such as "NDSolve" in Mathcmatica), and find the value of A that gives a periodic 
solution u, that is, that gives 

u(o) = u(2±*) = u(b) , it' (a) =u'(2±*) = u '(2±3fc) =u\b) =0, u (2±3*) = l, u{t) > Vt e [a, b] , 
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t 



t 



m = -2, A = m-0.8 m = -2, A = m - 1.0 m = -2, A = m - 1.2 




m = —3, A = m — 0.8 m = —3, A = m — 1.0 m = —3, A = m — 1.2 



FIGURE 3. The function u computed numerically via the first method. In each case, the 
graph is drawn for t G [a, with | = -a = 1.3108 (resp. 1.0, 0.8428) for m = -1 (resp. 
-2, -3). 



or equivalently, by the symmetries of V in Equation (15) and by the symmetries of u in Lemma 10, 
that simply gives 

(16) u'{a) = u'(2±*) = . 

Numerically finding this value of A for various values of to < 0, one finds that it is always to — 1, 
verifying the primary numerical result. 

Figure 3 demonstrates how this occurs, for the cases to = — 1, —2 and —3. In each case, the ODE 
solver produces a solution that satisfies Equation (16) precisely when A = m — 1. When A e (m — l,m] 
(for example, A = to — 0.8 as in Figure 3), the values of t (^ ^j^) closest to where u'(i) = 

satisfy |i - < ^p. When A e [to - 2, to - 1) (for example, A = m - 1.2 as in Figure 3), the 

values of t (^ ^±b) closest to ^2J± w herc u'(t) = satisfy \t - 22J±| > We need precisely 

u'i 2 ^) = u'(o) = 0, and this occurs exactly when A = to — 1. This has been checked for numerous 
other values of to < as well. 

8. Second method 

The second method we present here is more complicated than the first one, but we wish to consider 
it for the following two reasons: 1) it does not use an ODE solver and the algorithm is independent of 
the first method above, thus it gives a second fully-independent confirmation of the numerical result; 2) 
it is a stronger method in that it also gives estimates for other eigenvalues of £o, not just the first one. 

As we saw in Proposition 4 and the remark just after Proposition 6, in addition to Ao, —1 and are 
both also eigenvalues of £q for all to < 0, and —1 and are actually the second and third eigenvalues 
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of Co . The fact that this second method precisely estimates the second and third eigenvalues — 1 and 
gives us confidence that the method is accurately estimating Ao as well. 
We choose a basis {Bj — Bj(t)}°^ 1 for T as 

a 1 



\/b — a 



B i = \ll^-T.- cos [ 7:J l } t _ a a) ) l(,vov< ' 1J - ! ' :2Z ' 



This is an orthonormal basis for T with respect to the Euclidean L 2 norm on the interval [a, b] . Any 
u G T can be expanded as 

oo 

U = H a 3 B 3 
3 = 1 

for constants aj G ii. The Rayleigh quotient characterization (14) then gives 

, J2j.k>i aja k a jk 
A = mm — , 

Ei>i»?>o Ej,k>i°j a fc 



where 



n,, / -Bj-jr T Bi,.iU - I VB 3 B k dt. 



3 dt 

The integrals f°-BjjpB k dt can be explicitly computed as 



Bjj^B k dt - 



2 



4tt 2 



(b-af 



where Sj k is the Kronecker-delta function and [|] is the greatest integer less than or equal to |. 

Many integrals / VBjB k dt must be computed numerically (by a numerical integrator, such as "rein- 
tegrate" in Mathematica) , but for each n G Z + (n > 3) more than half of the entries of the nxn matrix 
(J^ VBjB k dt)™ k=1 can be determined to be zero simply by using the symmetry properties of the func- 
tions V and Bj. For example, V(a + t) = V(a — t) and B<2.{a + t) = ft (a — t), but Bs(a + t) = —Bz{a — t), 
and so f VB2Bzdt must be zero. 

Let us list the eigenvalues of £o in increasing order as 

A < Ai = -1 < A 2 = < A 3 < A 4 < ... -> +oo . 

Each eigenvalue appears in this list the same number of times as the dimension of its eigenspace. As 
noted in Proposition 5, liirij_ KX3 Xj = +oo. 

Remark. In fact, all eigenvalues have multiplicity at most 2. No eigenspace can contain three indepen- 
dent eigenfunctions, as Co is a second-order linear ODE. 

For each n G Z + , the nxn matrix (ajfc)™ fc=1 is symmetric, so it also has real eigenvalues, which we 
list in increasing order as 

X n) < A< n) <X ( 2 n) <...< X ( S . 

(n) 

Rayleigh quotient characterizations for Xj and Xj prove that 

(17) A? +1) > A? +2) > A? +3) > ... > Xj . 

Thus the limit lim^oo A^ exists and is greater than or equal to Xj. Using that {Bj = Bj(t)}j2. 1 is an 
orthonormal basis of J 7 , and so has dense span in jF with respect to the L 2 norm, further arguments 
with Rayleigh quotient characterizations give that in fact the limit is exactly equal to Xj : 
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Theorem 11. ([22]) lim^^ Xy' = \ 3 . 

Theorem 11 is proven in [22]. The essential facts behind its proof are well established and can be 
found in many sources ([3], [5], [26] for example), but we reference [22] because the result is given there 
in a situation exactly analogous to the one here. This result is a simple variant of the basic theory in 
finite clement methods (see the introductory chapters of [25] and [4], for example), and it is essentially 
only a variant of the Ritz-Galerkin method (see [18], for example). 

This theorem now provides us with our second numerical method for estimating Ao (and any other 
eigenvalue of Co) by simply finding the smallest eigenvalue (and other eigenvalues) of the matrix 
(Q!jfc)"fe = i for sufficiently large n. By the inequalities (17), it is clear that the estimates will be from 
above. Numerical results are shown in Table 1, and they again confirm the primary numerical result in 
this paper. 

Remark. Symmetries of the first eigenfunction can allow us to remove some Bj , if we are only looking 
for the first eigenvalue Ao- If u is the eigenfunction for the first eigenvalue, then the symmetries of u in 
Lemma 10 imply that u — X^>i a j^j with cij — when j > 2 is not an integer multiple of four. Then 
we can more quickly estimate Ao by using subspaces of T spanned by only B\ and B^k for k € Z + . 
However, this shortcut will not work for estimating the other eigenvalues of Cq. 







Lemma 


Lemma 


First 


Second 


Second 


Second 


Second 


Second 




b 
3 


7's lower 


9's upper 


method's 


method's 


meth.'s 


meth.'s 


meth.'s 


meth.'s 


m 


= —a 


bound 


bound 


estimate 


estimate 


est. 


est. 


est. for 


est. for 






for Ao 


for Ao 


for Ao 


for Ao 


for Ai 


for A2 


A3, A4 


^5, ^6 


-1/4 


2.0137 


-2.25 


-1.0522 


-1.25 


-1.245 


-0.992 


0.00543 


1.44 


4.45 


-1/2 


1.656 


-2.5 


-1.3643 


-1.5 


-1.4994 


-0.9987 


0.00102 


2.279 


6.75 


-1 


1.3108 


-3 


-1.9137 


-2 


-1.999 


-0.9993 


0.00062 


3.86 


11.021 


-2 


1.0 


-4 


-2.9483 


-3 


-2.999 


-0.9932 


0.00615 


6.94 


19.26 


-3 


0.8428 


-5 


-3.964 


-4 


-3.999 


-0.999 


0.00069 


9.943 


27.3 


-10 


0.4849 


-12 


-10.988 


-11 


-10.999 


-0.9974 


0.00256 


30.99 


83.46 


-20 


0.34683 


-22 


-20.994 


-12 


-20.999 


-0.983 


0.0168 


61.056 


163.61 



Table 1. Estimates for A and other eigenvalues of £ (computed using the values 
H = 1 and n = 13). 
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